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We investigate manipulation with Andreev bound states in Josephson quantum point contacts 
with magnetic scattering. Rabi oscillations in the two-level Andreev subsystems are excited by 
resonant driving the direction of magnetic moment of the scatterer, and by modulating the su- 
perconducting phase difference across the contact. The Andreev level dynamics is manifested by 
temporal oscillation of the Josephson current, accompanied, in the case of magnetic manipulation, 
also by oscillation of the Andreev states spin polarization. The interlevel transitions obey a selec- 
tion rule that forbids manipulations in a certain region of external parameters, and results from 
specific properties of Andreev bound states in magnetic contacts: 47r-periodicity with respect to the 
superconducting phase, and strong spontaneous spin polarization. 
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I. INTRODUCTION 

Recent advances in the development and experimental 
investigation of nanowire based Josephson junctions^'^'^ 
attract new attention to rich physics of mesoscopic 
Josephson effect. One of particularly interesting ques- 
tions concerns the possibility to employ Josephson quan- 
tum point contacts for quantum information processing. 
Such contacts contain a small number of generic two-level 
systems - Andreev bound levels, whose quantum states 
can be selectively manipulated and measured4iS By mod- 
ulating the phase difference across the junction one is able 
to induce the Rabi oscillation in the Andreev two-level 
system, and therefore to prepare arbitrary superposition 
of the Andreev states. A measurement of induced oscil- 
lation of the Josephson current allows for the Andreev 
level readout. Thus the pair of Andreev bound levels 
belonging to the same conducting mode may serve as a 
quantum hit^^ 

An interesting possibility to involve a spin degree 
of freedom in the contact quantum dynamics, and 
to use it for qubit application has been investigated 
by Chtchelkatchev and Nazarov.^ They considered a 
Josephson quantum point contact with spin-orbit inter- 
action and showed how to manipulate with the spin of 
the Andreev state. The properties of Andreev bound 
states in spin-active mesoscopic junctions and equilib- 
rium Josephson effect have been extensively studied in 
recent literaturej^ i-'^^'^^i^^'^'^i-'^'^'^^i^^i"'^^ non-stationary as- 
pects of interaction of individual magnetic scatterers with 
Josephson current have also been discussed i^^'^^i^°i^^ 

In this paper, we investigate the methods of manip- 
ulation with the Andreev states in Josephson quantum 
point contacts containing a magnetic scatterer, e.g. mag- 
netic nanoparticle situating between the superconducting 
electrodes.-^^^^^ We investigate two manipulation meth- 
ods: (i) time variation of the superconducting phase 
across the contact, and (ii) time variation of the direc- 



tion of magnetic moment of the scatterer. We find that 
in both the cases the Josephson current exhibits Rabi 
oscillation under the resonant drive within a certain in- 
terval of biasing superconducting phase. In the case of 
magnetic manipulation, the effect may only exist if the 
Andreev states are initially spin polarized, then the cur- 
rent oscillation is accompanied by oscillation of the An- 
dreev states spin polarization. The phase interval, where 
the Rabi oscillation can be excited, decreases with in- 
creasing strength of the magnetic scatterer, and eventu- 
ally disappears at large enough strength; in particular in 
the TT-j unction regime, the Rabi oscillation is completely 
forbidden. This selection rule results from specific prop- 
erties of the bound Andreev states in magnetic junctions 
as we will show. 




FIG. 1: Sketch of a magnetic Josephson point contact: su- 
perconducting reservoirs are connected by a nanowire of 
the length L smaller than the coherence length, magnetic 
nanoparticle creates local classical magnetic field. 



For a static scatter the spin rotation symmetry around 
the direction of its magnetic moment is preserved. This 
allows for the contact description in terms of a two- 
component Nambu spinor,^* similar to non-magnetic 
junctions, thus avoiding a double counting problem. 
Within such an approach, the two bound Andreev lev- 
els per conducting modes are only relevant, giving com- 
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plete quantitative description of the stationary Joseph- 
son effect as well as the resonant two-level transitions 
and non-stationary current response under the phase 
manipulation.- Consideration of the spin conjugated 
Nambu spinor gives a completely equivalent physical de- 
scription in terms of a reciprocal pair of Andrccv bound 
states; both the pictures mirror each other. 

Time variation of the direction of the magnetic mo- 
ment of the scatterer leads to a violation of the spin 
rotation symmetry, and induces coupling between the 
spin conjugated Nambu spinors. This results in a uni- 
tary rotation in the extended space of the four Andreev 
bound states. It turns out, however, that this rotation 
splits into two equivalent rotations in invariant two-level 
subspaces, which mirror each other. Thus the contact 
response can also in this case be explained in terms of 
the two- level Rabi dynamics. For the two different ways 
of magnetic manipulation considered - instant switching, 
and small-amplitude resonant oscillation of the direction 
of the magnetic moment of the scatterer, the Andreev 
two-level dynamics has a physical meaning of precession, 
and nutation of the spin polarization of the Andreev lev- 
els, respectively. Thus the spin polarization of Andreev 
states is required in order to observe a non-trivial dynam- 
ical response.— Such a possibility naturally exists, as we 
will show, in the contact under consideration: the equi- 
librium Andreev states exhibit strong spin polarization, 
up to the maximum values ±1/2 at low temperature, 
in certain regions of the superconducting and Zeeman 
phases. 



II. CONTACT DESCRIPTION 

Consider one-mode quantum point contact with super- 
conducting electrodes connected by a normally conduct- 
ing nanowire as shown in Fig. [T] The left and right 
electrodes (L,R) are described with the BCS Hamilto- 
nian, 



Hs ^ [ dx y2 iplix) Kx,t)] tpaix) 

+A*(a;, t)i^^ix)i^i{x) + A{x, 



Aiming to investigate the effect of time variation of 
the superconducting phase difference ip across the junc- 
tion, we consider the order parameter having the form, 
A{x,t) = Ae"^(^^*)/2^ ip{x,t) = ip{t) sgnx, and the elec- 
trochemical potential, fi{x,t), having the form, ii{x,t) — 
Ep — which provides the electro-neutrality con- 

dition within the electrodesj^- 

To explore the spin properties of the Andreev states, 
we assume that the contact nanowire contains a mag- 
netic scatterer, e.g., magnetic nanoparticle (Fig. [1]). We 
assume for simplicity that the magnetic field H(a;, t) in- 
duced by the scatterer is localized within the nanowire on 
a distance I smaller than the distance L between the elec- 
trodes, thus not affecting the superconductivity within 



the electrodes; furthermore, we will treat this magnetic 
field as a given external parameter neglecting the back 
action effect from the current. Then the Hamiltonian of 
the normal region of the junction has the form. 
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where U{x) is the scalar potential of the scatterer. We 
assume a symmetric with respect to a; = spatial dis- 
tributions of the scalar potential and the magnetic field, 
and also a fixed direction of the magnetic field, which will 
vary with time during the manipulation. 

In the stationary case, the Hamiltonian ^ preserves 
the spin rotational symmetry around direction of the 
magnetic field. Choosing the spin quantization axis 
along this direction, we describe the electron propaga- 
tion through the normal region of the junction with a 
transfer matrix, To, 



.(/5/2) 



(3) 



where d — diag(d|,dj^) and f — diag(r|,r|). Contact de- 
scription in terms of the spin active scattering matrix has 
been extensively discussed in literature f^i2^ The impu- 
rity scalar potential produces spatially symmetric scat- 
tering with transmission amplitudes, and reflec- 
tion amplitudes, ''|,r|, which may be different for dif- 
ferent spin orientations (spin selection). The scattering 
phase shift [3 between the opposite spin orientations is 
induced by the Zeeman effect. 



/3 = 



(4) 



The physical observables of interest, Josephson cur- 
rent and spin polarization, are described with a single 
electron density matrix p{x,x',t) associated with a two- 
component Nambu field^^ \I'(x,t), 

^{x, t) = jl^^ JJ^ , p(x, x\ t) ^ {^{x, t)^\x', t)) 

(5) 

(here the angle brackets indicate statistical averaging). 
Alternative description is given by a density matrix, 
p{x,x',t) associated with the spin conjugated Nambu 
field, *(a;,i), 

§(x,t) = - ~p{x,x',t) = {^{x,t)¥{x',t)). 

(6) 

The two Nambu fields are connected via a fundamental 
symmetry relation imposed by the singlet nature of the 



BCS pairing, 

^{X,t)^l<jy{¥f . (7) 

In what follows we will use these two reciprocal rep- 
resentations, (/(-representation and ^-representation, re- 
spectively (see Appendix B). This will allow us to avoid 
the redundancy, introduced by Eq. ([7]), of a commonly 
used four-component Nambu formalism )^'^°i^^i^^i^^'^"^i^^ 
and to explicitly show that the stationary Josephson ef- 
fect as well as non-stationary response to the phase ma- 
nipulation can be fully understood within the framework 
of the two bound Andreev states per conducting mode 
associated with either of these two reciprocal representa- 
tions. 



III. ANDREEV STATE SUBSYSTEM 




FIG. 2: Andreev level energy spectrum for sin(/3/2) < ^/D 
{D=0.9, /3=0.5), upper panel; 9- functions, shown on the lower 
panel, define the spectrum discontinuity points at 2nn± ipo- 
Vertical dotted lines separate regions of weak and strong Zee- 
man effect. 



In this section we outline the properties of the sta- 
tionary Andreev states, which are of importance for the 
discussion of non-stationary effects. The details of deriva- 
tions are presented in Appendix A. 

Within the two-component Nambu formalismf2^ the 
quasiparticle states in the stationary contact are given 
by the Bogoliubov-deGennes equation. 



(8) 



supplemented with the boundary condition at the con- 
tact, which for a short contact is given by the transfer 
matrix, T = ex.\){i(j zip / 2)Te, Eq. (|A3p . The energy spec- 
trum of the bound states consists of the two levels, ac- 
cording to Eq. (jA9p . 



9sAcos{sr]- P/2), 



(9) 



where Og — sgn[sin(s77 — /3/2)], and parameter rj is de- 
fined through equation, sin?7 = \/15sm((p/2), D is the 
contact transparency. This energy spectrum is asymmet- 
ric with respect to the chemical potential, E = 0, and is 
47r-periodic with respect to the phase difference as shown 
on Fig. [21 the spectral branches cross at (p — 27rn. On 
the figure are also shown the energy levels of a recipro- 
cal Andreev level pair, Es = —Eg, whose wave functions, 
4>s = {i<^y)4>*s^ ^'^^ associated with the spin conjugated 
Nambu field, Eq. ^ (^-representation, see Appendix B). 



The reason for the 47r-periodicity and the level cross- 
ings is the symmetry of the problem which is reflected in 
the property, Eq. (jA12p . Introducing the parity operator 
P that permutes the wave functions at the left and right 
sides of the junction, P(t){x < 0) = (t>{x > 0), we write 
Eq. (|A12P on the form^ 



,4>six), A = P 



CT. 



(10) 



According to this equation, the Andreev level wave func- 
tions are simultaneously the eigen functions of the sym- 
metry operator A with eigenvalues sOg- Furthermore, 
the energy branches, E- and Ej^, in the neighboring 
phase intervals in Fig. [2l correspond to different eigen 
values of the operator A. The properties of the An- 
dreev states are therefore qualitatively different within 
the phase intervals, where dg — ^'-s, and Og = —O^g. 
To emphasize the difference we will refer to the for- 
mer ones as the regions of strong Zeeman effect (ZE), 
sin/3/2 > \fDsmip/2, and the latter ones as the regions 
of weak ZE, sin/3/2 < ^/D shnp/2. These regions are sep- 
arated by the points where the energy levels touch the 
continuum, ip — 27m ± Lp^, po = 2arcsin(l/-\/Dsin/3/2). 
At sufficiently large Zeeman phase, sin/3/2 > ^/D the 
weak ZE regions disappear, while strong ZE regions 
spread over the whole superconducting phase axis. The 
bound energy levels, E±{p), depart from the continuum 
forming " cigars" , see Fig. [3] , they belong to the orthogo- 
nal eigen subspaces of the operator A at all phases. The 
TT-contact is realized in this regime, at /3 = tt^^ when the 
reciprocal cigars coincide and situate symmetrically with 
respect to E — 0. 



A qualitative difference between the regions of strong 
and weak ZE is illustrated by the properties of Josephson 
current. Starting with a general expression for the charge 
current through the density matrix, 

(11) 

we truncate it to the Andreev level subspace, 

lA{t) = J2 Iss' (I 5g,g ~ pg'g{t)\ . (12) 
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together with the spectrum discontinuities result in a pe- 
culiar phase dependence of the spin polarization of An- 
dreev levels. The z-component of electronic spin density 
in the contact is given by equation, 



x') — Trp{x^ x' , t) ] 



(16) 



Truncating this equation to the Andreev level subspace, 
and integrating over x using normalization condition for 
the bound state wave functions, we find the spin polar- 
ization of the Andreev level pair, 



5a = (1/2)(1 -/+-/_), 



(17) 



FIG. 3: Andreev level energy spectrum for sin(/3/2) > ^/D 
{D—0.1, P = 2.8). Strong ZE region spreads over whole phase 
axis. 



Here 



dxdx'Ti ^p{x, x' ,t)(j)s' {x')(f>l{x)'^ 

(13) 



is the density matrix in the Andreev level representation, 
Eqs. (jAlip . (jA12p . the trace refers to the electron-hole 
space. The current matrix Iss' reads. 



_2e 
n 







d^E- 



+ {1~ 9s0-s)VRDsm 



if VOC^ fo 1 



2e 



1 



(14) 



The diagonal elements give the expectation values of 
the currents of Andreev states, while the off-diagonal 
part describes the current quantum fluctuation.'' The off- 
diagonal part is finite in the weak ZE regions, where 
9s9-s — —1, while it is zero in the strong ZE regions. 
This implies that the current quantum fluctuation is fully 
suppressed in the strong ZE regions: here the current 
matrix commutes with the Andreev level Hamiltonian, 
which is diagonal in this representation hgs' — Sss'Es- 

In the equilibrium, pss' = 5ss'nF{~Es), and the An- 
dreev current in Eq. (|12p takes the form, 



where fs — Pss are the population numbers. 

Thus the spin polarization of the Andreev levels is en- 
tirely determined by their (generally non-equilibrium) 
population numbers: For empty Andreev level pair, 
/± = 0, the spin polarization is 5 = 1/2, while for fully 
populated levels, /± = 1, it is 5 = —1/2. For the sin- 
gle particle occupation of the level pair, /+ + /- = 1, 
the spin polarization is zero, 5 = 0. In non-magnetic 
contacts the equilibrium spin polarization of Andreev 
levels is always zero, S" = 0, by virtue of the identity, 
ni?(— i?_|_) -I- npi—E^ — 1, that holds due to the spec- 
trum symmetry, E+ = —E^. In magnetic contacts, the 
spin polarization sharply varies in the f — P parameter 
plane, see Fig. 21 At zero temperature, fs — 0{Es), the 
spin polarization is zero in the regions bound by the lines, 
sm{(p/2) = sm{(3/2)/VD and sin((/j/2) = cos(/3/2)//D; 
the first region corresponds to a weak ZE. Outside these 
regions the Andreev levels are strongly polarized, S = 
1/2. At small contact transparencies, D < 1/2, the above 
mentioned lines do not overlap, while at D > 1/2 they 
do, forming the island of strong negative polarization, 
S = —1/2, around the point, (p = tt, ^ = 7r/2, see Fig. S) 
This island grows with increasing transparency eventu- 
ally touching the lines, /3 = and j3 = it a,t D — 1. 



Ia 



eA 



Y^d^Es tanh(S,/2r), 



(15) 



One can confirm by direct calculation that the current 
of the continuum states vanishes, thus Eq. (|15|1 repre- 
sents the total Josephson current. This equation coin- 
cides with the result obtained using the four-component 
Nambu formalism,— li^iiiii^ii^ii^iifi and it can be also ob- 
tained by working with the reciprocal Nambu represen- 
tation using equation (|B9p . 

Spin polarization of the Andreev states plays an im- 
portant role in magnetic manipulation. Asymmetry of 
the Andreev spectrum with respect to the zero energy 




FIG. 4: Equilibrium spin polarization of Andreev levels at 
zero temperature and at D > 1/2; the central region with 
S — —1/2 emerges at D = 1/2 a.t ip — tt, l3 = n/2, and grows 
with D reaching the lines /3 — and P = n at D — 1. 



5 



In contrast to the charge density, the spin density in 
Eq. obeys the conservation equation, dtS{x,t) + 

dxls{x,t) — 0. The consequence of this is the zero spin 
current of the Andreev states, Isa — 0: under the sta- 
tionary condition the partial spin current of the Andreev 
state must be constant in space, and being proportional 
to the bound wave function it vanishes at the infinity; 
thus it is identically equal to zeroi^ 



IV. PHASE MANIPULATION 

Now we turn to discussion of the Andreev level dy- 
namics under the time dependent phase. Similar to non- 
magnetic junctions this dynamics involve only the bound 
levels belonging to the same Nambu representation, "^^i^ 
as shown on Fig. [S] The frequency of the phase time 
variation must be small compared to the distance to 
the gap edges to prevent the level-continuum transitions, 

Lu < A - |£;,|. 

Time evolution of the contact density matrix, p{t), is 
governed by the Liouville equation, ihdtp = [h, p], with 
the Hamiltonian of Eq. ([8]), supplemented with the non- 
stationary boundary condition, T = e'xjp(iaz^{t)/2)Te. 
Now we truncate the density matrix using the instanta- 
neous Andreev eigenfunctions, (jisifit)] (cf. Eq. p^ ). 



Pss'it) = (MtMtMs'it)). 



(18) 



This density matrix obeys the Liouville equation ihdtp — 
[H, p] with a truncated Hamiltonian, which in this basis 
is given by equation, 

i/ss' = {(t>s\h\(j)s')-{(t>s\'itidt\4>s') = Es5ss'~fiv{'t>s\'id^\4>s'), 

(19) 

where Es[ip{t)] and 4)s[f{t)] are given by Eqs. (|A9p and 




FIG. 5: Interlevel transitions induced by time oscillation of 
the phase difference j_ shadow regions indicate forbidden re- 
gions; transitions in <j) representation (reduced- intensity lines) 
are equivalent to the transitions in (^-representation (full- 
intensity lines); 73=0.9, /3=0.5. Inset: transition matrix el- 
ement as function of ip. 



The matrix element, {(t)s\id^p\(t)s') ^ is found to be zero 
for s' = s, while for the interlevel transitions, s' = — s, it 
reads, 



A(vp) 



= is{l - dsO^s) 



AM 
2 ' 



RD 



■ 2 f VC+C- 

sm — - 



(20) 



2e2(C++C- 



From this we find that in the weak ZE regions, OgO-s = 
— 1, the non-stationary Andreev level Hamiltonian has 
the form, 



Hit) 



E+ 








fi'PH'p) , /D sin I > sin ^. 

(21) 

This equation provides generalization to a magnetic junc- 
tion of the Hamiltonian derived in Ref. @ : the level cou- 
pling coincides with the one in non-magnetic junction 
when /3 = 0, and remains finite when p ^ but only 
inside the weak ZE region, decreasing towards the edges 
of this region, see inset in Fig. \5\ In the strong ZE region 
the matrix element is identically zero [OgO^s = +1)) find 
the Hamiltonian is diagonal. 



H{t) 



E+ 






E_ 



. .13 
D sin — < sin — . 



(22) 



Thus we conclude that no operation with Andreev levels 
is possible in the strong ZE regime. 

Eq. ([2T|l is convenient for calculation of Rabi oscillation 
in the Andreev level system under the resonant driving, 
</?(t) = Lp -\- Jsinwt, Lo = (i?_|_ — E^)/h. Inserting this in 
Eq. jUl), we get 



H = — a, 
2 



HASlU cos Ujt (T,, 



(23) 



Assuming small amplitude of the phase oscillation, S <C 
1, and using rotating wave approximation, we find the 
time dependent density matrix in the rotating frame, 

p{t) = p{0) - M^^-f-i^) _ cosASujt ~ sin ASLuta^), 
^ (24) 

where p(0) = diag{f+{0), f-{0)). Rabi oscillation of the 
Andreev levels generate a time dependent Josephson cur- 
rent, 

m = m - I (/+(0) - /-(0))sin2 ^ Y.9<fEs . 

(25) 

The Rabi oscillation, Eq. (|24|) . and the time oscillation 
of the Josephson current, Eq. ([^5]) . vanish if the Andreev 
levels are initially fully spin polarized, S{0) = ±1/2, since 
/+(0) = /-(O) in this case. Thus an additional require- 
ment for the phase manipulation is to bias the contact in 
the region outside the negative-spin island on Fig. U) 
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V. SPIN MANIPULATION 

Now we proceed with the discussion of the spin ma- 
nipulation. We consider the two ways of driving An- 
dreev level spin, presented in Fig. [S| (i) rapid change of 
the direction of the magnetic moment of the scatterer (dc 
pulsing), and (ii) harmonic oscillation with resonance fre- 
quency of the magnetic moment direction (rf-pulsing) . In 
both cases the spin rotation symmetry is violated, there- 
fore the junction dynamics cannot be described with only 
one Nambu pseudo-spinor, but involves both the spin 
conjugated Nambu pseudo-spinors. The interlevel tran- 
sitions in this case physically describe a rotation of the 
Andreev level spin. 




e(t)«i 

FIG. 6: Sketch of manipulation with magnetic field; a) instant 
switching of direction of magnetic field; b) small-amplitude 
oscillation of magnetic field with resonant frequency (electron 
spin resonance) 



A. dc pulsing 

Suppose the Andreev levels are initially prepared in a 
stationary state with non-zero spin, which points along 
the applied magnetic field (z-axis) . Such states were dis- 
cussed in previous sections. Let us now suppose that the 
magnetic field is rapidly rotated by angle around y-axis, 
as shown on Fig. [6^. Such a manipulation is described 
by rotation of the electronic T-matrix, in Eq. ([5]), 

Te UT,U\ C/ = f i '^"^ i V (26) 
V-sm| cos|; 

and it mixes the Nambu pseudo-spinors 5" and ^. To 
describe the effect of this manipulation, we introduce the 
extended four-component Nambu space, (^',4')"^, and 
the corresponding single particle density matrix, 

(27) 



This density matrix operates in the Hilbert space 
spanned by the extended eigen basis, 

$.(.) = (^f)) , M^) = (^^^,).) . (28) 

The transformation U induces rotation of the extended 
basis 

^a{x) ^ U^a{x), ae{iy,i)}. (29) 

The eigen energies, however, remain the same since such 
a rotation just corresponds to changing the spin quanti- 
zation axis. 

The wave functions, Eq. ((28l) . form a complete set in 
the extended space, and thus any operator can be ex- 
pressed through them, 

Aix, x') = J2 <^>c.ix)^lix')A^p. (30) 

In particular, for a stationary system with spin rotational 
invariance we have for the density matrix, 

n(a;,a;') = ^<i>a(a;)<i>J.(x')/a, (31) 

a 

and the Hamiltonian, 

ff(x) =^$„(a;)$t(a;)i;^. (32) 

a 

However, one has to remember that this description is re- 
dundant, and rigorous constraints hold on the occupation 
numbers, f^ = l — f^, and eigen energies, — —E^. 

We now write down the Hamiltonian after the mag- 
netic field rotation on the form, 

H{x) = U'^c.{x)E^<^l{x)U\ t > 0. (33) 

a 

In the initial basis, this Hamiltonian is represented with 
the matrix, 

or explicitly, 

H^^, = cos^ -E^6^^, - sin^ - ^ {(l)^\(f>^,)E^{4>^\(p^,) 
Ho9- = - cos^ -^EJ^^, + sin^ - ^ {<PA<^^^^)E^,{(t)^,\(i)u') 
Hi,i^: = -sin6'((/),,|(^,y>^^^^-t-^. 

(35) 

Here the orthogonality relations were used, (^i/j^p) = 
(<^^|0^) = 6^^, and (<?!)^|0^) = 0, Eq. 
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At this point we restrict ourselves to the Andreev level 
subspace, and present the truncated Hamiltonian on the 
form (using the symmetries Eg — —Eg and {(f>+\(f)^) — 



we obtain, 



ij(4) = 

where 

Eo,W 
V 



/Eo + W 


V 











V 

~Eo + W -V 
-V -Eo -W 
Eo-W/ 



cos^ - -sin^ -iMp 



-Msm9 



E. 



E_ 



(36) 



(37) 



and the inter level matrix element, 

M=(0+|0_)-((?, -0_,)CO! 



VC^- 

2 C+ + C- 



(38) 



The matrix element psp equals to zero in the strong 
ZE region {9s — 0-s), thus the manipulation does not 
produce any effect there, which is similar to the phase 
manipulation, see inset in Fig. [T] 

The Hamiltonian ([36|) has a block-diagonal form de- 
scribing identical rotations in the two orthogonal sub- 
spaces, spanned with the eigen vectors (</)+, </)-), and 
(0_, ^4-). Thus the problem reduces to solving for two 
physically equivalent two-level systems. Choosing the 
subspace {4>+, </>-), we have the two- level Hamiltonian, 



W V 



V -wj ■ ^^^^ 

Introducing projection operators on the eigen subspaces. 



A=± 



= + f cos2 ^ + sin2 ^|Mp 



(40) 



where XhQ are the eigen energies, we have for the time 
evolution of the two-level density matrix. 



V 

a(t) = —i — sin Jit, b{t) = a(t) 
nil 



W 

cos ilt + i — sin fit 



(43) 



where 5a(0) = (l/2)[l-/+(0)-/_(0)] is the initial spin- 
polarization of the Andreev levels as given by Eq. (|17p . 
Thus no rotation is induced for spin unpolarized An- 
dreev levels. Furthermore, the frequency of the rotation 
is proportional to the level splitting, hfl cx E^ — E^ = 
E^ + E-, i.e. to the magnetic field. 

The time-evolution of the occupation numbers fs (t) of 
the Andreev levels in the </)-representation is extracted 
from Eq. (gS]), 

Mt) - /,(0) + 2^^(0)^^ sin2 nt. (44) 

This relation illustrates the non-unitary evolution of the 
Andreev levels in this representation, f+{t) + f-{t) ^ 
const. 

Eq. (UH) allows us to obtain the time dependence of 
the spin polarization of the Andreev levels. 



SA{t) = l-{l-U{t)-f-{t)) ^ Sa{0) [I 



2V' 

Jmj' 



■ sin Qt 



(45) 



To calculate the Josephson current, we use the expres- 
sion through the current matrix in the 0-representation, 

Eqs. (na), (m. 



(46) 



The diagonal elements of the density matrix, pss, are 
given by Eq. (UH). On the other hand, the off-diagonal 
elements, pss' , equal zero because the spin manipulation 
does not induce transitions between the eigen states of 
the same (either (jj- or (jj-) representation. Therefore, 



2e 



S 

= /(O) + ^5^(0) ^ sin^ nt J2 d^Es 



(47) 



n(2)(t) e-*'^n(2)(0)e^^ = ^e'(^'-^)^"PAn(2)(0)PA'. 

AA' 

(41) 

Assuming the initial density matrix to be stationary (not 
necessarily equilibrium), and expressing it through the 
level occupation numbers of the ^-representation, 

""'^^(TlJ^.O)). (42) 



Summarizing, the conditions for the observation of the 
non-stationary contact response, biasing in a weak ZE re- 
gion with a finite spin polarization, can be only fulfilled in 
the central island region on Fig. [4] with negative polariza- 
tion. This constraint can be relaxed by pumping initial 
level populations away from equilibrium as suggested in 
Ref.8. 

Concluding this section we note that the dc pulsing 
of the magnetic field does not allow one to reach every 
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E 




FIG. 7: Interlevel transitions induced by magnetic manip- 
ulation; shadow regions indicate forbidden regions; transi- 
tions between the levels {E-, _E+) (reduced- intensity lines) are 
equivalent to the transitions {E+,E-) (full- intensity lines); 
D—0.9, p—0.5. Inset: transition matrix element as function 
of ip. 



the density matrix in terms of these instantaneous eigen- 
functions, 

nix,x',t) = Y,umo.ix)<Pl{x')uHt)n^fj{t). (50) 

The matrix Ilap{t) satisfies the Liouville equation with 
the Hamiltonian, 



Hc^pit) = J dx<^i{x)UHt) [H{x) - ihdt] U^fiix) 
=E^&c.p - ih{<^^\U\t)dtU(t)\<^p). 



Inserting Eq. we get, 



ihu\t)dtU{t) = -htO ' 



2' \i 



(51) 



(52) 



Truncating to the Andreev level sub-space we have the 
Hamiltonian, 



point on the Bloch sphere: this is due to the fact that U is 
not an invariant operation on the Andreev level subspace 
of the extended Nambu space, which physically means a 
leakage to the continuum (similar effect exists also for 
the phase manipulation with dc pulses^). However, for 
small rotation angles and not close to the edges of the 
weak ZE region, the matrix element M is close to unity, 
and the leakage is small; it can be further reduced by 
using rapid adiabatic change of the magnetic field, i.e. 
rapid on the time scale of the Andreev level splitting but 
slow on the time scale of the distance of the Andreev 
levels to the continuum. This shortcoming does not exist 
for the resonant rf-pulsing. 



B. rf-pulsing 

Now let us consider a time-dependent rotation of the 
magnetic field, 

T ^u{e{t))Tu\e{t)). 

where 0(t), as before, is the angle of rotation around the 
y-axis, see Fig. |6b. We can now define the instantaneous 
eigenstates 



(48) 



satisfying the instantaneous boundary condition, Eq. ([3]), 
with Te UTgW . The time-dependent Hamiltonian can 
then be written similarly to Eq. p3|) . 



Hix, t)=J2 Um^{x)E^<i>lix)UHt). (49) 

a 

Since the energy eigenvalues do not depend on the 
direction of the quantization axis they remain time- 
independent. Now similarly to Eq. pip we can expand 



H = 




, g^-dtOM. (53) 



This Hamiltonian can again be presented in a block- 
diagonal form describing two equivalent two-level sys- 
tems. Choosing the subspace spanned by ((/)+,0_) and 
driving the magnetic field at exact resonance, uj = — 
E^ = E^ + with small amplitude, d{t) = Oosinuit, 
6*0 <C I, we have the two- level Hamiltonian, 



E^ —iM^r cos ujt\ n —^Q Af 

ihnrCOSLOt -E_ J' ^ir-^Oo^M, 

(54) 

from which we obtain the Rabi oscillation of the popula- 
tion numbers of the (/)-representation, 

Uit) = /+(0) cos^ ^ + [I - !-m sin^ 



[I - !.{t)\ = [I - /_(0)] cos^ ^ + UiO) sin^ 



2 



(55) 



or introducing explicitly the the Andreev level spin, 
/,(t) = /.(0) + 25^(0)sin2^. 



(56) 



This equation again illustrates the non-unitary evolution 
of the Andreev levels in the (/^-representation. The time 
evolution of the spin polarization, and the Josephson cur- 
rent then become, respectively, 

m =/(0) + ^^5^(0) sin^ 'Jf Yl 
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VI. DISCUSSION 

Manipulations with the Andreev levels generate 
strongly non-equilibrium states whose lifetime is re- 
stricted by relaxation processes. Let us qualitatively 
discuss the relaxation mechanisms relevant for the non- 
equilibrium states induced by manipulation methods dis- 
cussed. 

The phase manipulation affects the difference of the 
populations of the Andreev levels belonging to the same 
Nambu representation while keeping the total population 
the Andreev level pair unchanged. At zero temperature, 
and for relatively small frequency of the qubit rotation 
compared to the superconducting gap, the states of the 
continuum spectrum are either empty or fully occupied, 
and therefore the exchange between the continuum and 
the Andreev levels is exponentially weak.— Therefore the 
relaxation predominantly occurs within the Andreev level 
system. In the strong ZE regions, the interlevel relax- 
ation caused by interaction with electromagnetic envi- 
ronment should be suppressed due to the vanishing tran- 
sition matrix element, Eq. (j20p . i.e. for the same reason 
that prevents the phase manipulation. One may expect 
to prolong lifetime of the excited states by taking ad- 
vantage of this property and shifting adiabatically the 
phase bias into the strong ZE region after the manipu- 
lation has been performed. Such an operation, however, 
requires a passage through the one of the singular points, 
tf = 27rri ± (y3o, where the Andreev levels touch the con- 
tinuum; at this points the quantum state escapes in the 
continuum and quantum information is lost. 

The absence of the interlevel relaxation in the strong 
ZE regions has interesting implications for the obser- 
vation of 47r-periodicity of the Andreev level spectrum 
discussed in Section III. The equilibrium Josephson cur- 
rent, Eq. p5|) . is 27r-periodic and does not reveal the 
47r-periodicity property of the Andreev states. This may 
change under non-equilibrium condition, when a small 
voltage is applied to the junction. In this case, supercon- 
ducting phase becomes time dependent, ip = 2eV/h, and 
the Andreev levels adiabatically move along the i^-axis 
keeping constant level population during long time (lim- 
ited by a weak level-continuum quasiparticle exchange). 
If the magnetic effect is weak while contact is trans- 
parent, sin(/?/2) < a/D, the levels touch the contin- 
uum every Josephson cycle, and the level population will 
be periodically reset^ leading to the 27r-periodicity of 
the Josephson current. However, in the "cigars" regime, 
sin(/3/2) > -s/D, depicted on Fig. [3l the levels are iso- 
lated from the continuum, and the level population may 
remain unchanged during the time greatly exceeding the 
Josephson period. This will lead to the the 47r-periodicity 
of the ac Josephson current, and could be experimentally 
detected by observing anomalous Shapiro effect with only 
even Shapiro steps present. The effect would be the most 
pronounced for the 7r-junctions, /3 = tt. Similar effect has 
been discussed in a different context of unconventional 
superconductor junctionsi^ 



Manipulation with the Andreev level spin affects the 
spin polarization of the Andreev levels, and thus, at first 
glance, a relevant relaxation mechanism would require 
some spin active scattering. Since magnetic interactions 
in superconductors are usually rather small compared to 
non-magnetic interactions, e.g. with electromagnetic en- 
vironment, one would expect a long lifetime of Andreev 
spin excitations.^ However, one should take into account 
the relation between the spin polarization and population 
of the Andreev level pair Eq. (|B9[) : The non-equilibrium 
spin polarization is associated with non-equilibrium pop- 
ulation of the Andreev level pair, which can be relaxed by 
any non-magnetic interaction. Consider, for example, the 
process of approaching the equilibrium state in magnetic 
contact after the phase bias has been suddenly changed. 
This will first create a non-equilibrium state in the An- 
dreev level system, both in terms of individual level popu- 
lations, and total population of the level pair, which then 
will rapidly relax to local equilibrium via the interlevel, 
and level-continuum quasiparticle transitions induced by 
(presumably) strong non-magnetic interactionii Such in- 
teraction does not change the total spin polarization of 
contact electrons, since it preserves the spin rotation 
symmetry, however it is able to transfer the polarization 
from the Andreev levels to the continuum states. The 
total polarization is maintained in the local equilibrium 
by shifting the energy argument in the Fermi distribu- 
tion function by an energy independent constant. This 
constant will slowly relax to the zero value in a second 
relaxation stage, due to spin-flip processes. 

Thus we conclude that decoherence of the states gen- 
erated by the magnetic manipulations basically results 
from the same physical interactions that destroy excited 
states produced by the phase manipulations, and there- 
fore one should not expect significant differences of the 
respective lifetimes. 



VII. CONCLUDING REMARKS 

In conclusion, we studied the properties of Andreev 
bound level system, and various ways of manipulation 
with them in Josephson quantum point contacts contain- 
ing magnetic scatterers. 

In practice, such contacts can be realized by attach- 
ing magnetic nanoparticles or molecules to the contact 
bridge;-^ another possibility is to insert magnetic macro- 
molecules, e.g. doped metallofuUerene in the contacti^ 
Coulomb blockade regime in molecular dots offers ad- 
ditional possibility due to uncompensated spin of odd 
electronic configurations on the dotFii^i^ii^ 

In the studied cases of resonant driving the super- 
conducting phase difference and the direction of mag- 
netic scatterer, the contact response consists of a time- 
oscillation of the Josephson current, and for the magnetic 
drive, also oscillation of the Andreev level spin polariza- 
tion. We identified the regions of external parameters 
where these oscillations can be excited. The correspond- 
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ing selection rule results from specific symmetry proper- 
ties of the bound Andreev states in magnetic contacts: 
47r-periodicity of the level spectrum and strong sponta- 
neous spin polarization. 

In all the studied cases, the non-stationary contact re- 
sponse results from resonant dynamics of two physically 
identical two- level systems (for one conducting mode), 
whose evolutions mirror each other. This is the man- 
ifestation in a non-stationary regime of the redundancy 
(double counting) of the four-component Nambu descrip- 
tion of the Josephson effect in magnetic contacts. Due 
to the fundamental constraint, Eq. ([7|), imposed by the 
singlet pairing the four-component Nambu field possesses 
the algebraic structure of Majorana fermion^^^ thus de- 
scribing only two physical degrees of freedom rather than 
four. In contacts with magnetic impurities under discus- 
sion, these two physical degrees of freedom relevant for 
the stationary Josephson effect correspond to the two 
Andreev bound states per conducting mode, similar to 
the case of non-magnetic contactsj^'^'^ 
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APPENDIX A: BOUND STATE WAVE 
FUNCTIONS 

To explicitly construct the wave functions of the An- 
dreev states, we consider a quasiclassical approximation 
for (j){x) by separating rapidly oscillating factors e**'^''^, 
and slowly varying envelopes (j)^{x), 

(j){x) ^ (/)+(x)e''=^^ + 0-(x)e~*'^'^^. (Al) 

The envelopes, 0^, satisfy a quasiclassical Bogoliubov-de 



Gennes (BdG) equation 



37 



(A2) 

Furthermore the superconducting phase can be elimi- 
nated from this equation and moved to the boundary 
condition by means of the gauge transformation, (j) — > 
exp ( i cr^ sgn X / 4 )(/) . 

The boundary condition at the contact for the qua- 
siclassical envelopes, 0^(0), follows from the electronic 
transfer matrix in Eq. ([3]). We assume for simplicity the 
short contact hmit, L <C Co, where is the supercon- 
ducting coherence length, thus neglecting the energy dis- 
persion of the scattering amplitudes. Then it is easy to 
establish that the transfer matrix for holes has the same 
form as for the electrons. Thus the boundary condition 



connecting the left (L) and right (R) electrode wave func- 
tions can be written on the form, 



<!>- 



(A3) 



Elementary solutions to a stationary BdG equation, 
{±vpp(7z + AcTj;) (j)^ = Ecjr^^ have the form for given en- 
ergy \E\ < A, 



1 / ±ia7/2\ 



where 



V2 ye 



E 



cos7^— , sm7 = -^, C = - £:2. 



(A4) 



(A5) 



Index a is defined by the zero boundary condition at the 
infinity. The matching condition, Eq. (|A3[) then reads. 



— g*o-2(v/2)2i 



B-(h- 



(A6) 



where the coefficients , are to be determined by 
this equation and the normalization condition. The solv- 
ability of this matching requires. 



cos(27 + 13) = R + Dcosip, 



(A7) 



where R = r-^r^ and D — d^d^ play the role of effective, 
spin-symmetric reflection and transmission coefficients^ 
(cf. Ref. where a more general form of this equation 
has been derived). Introducing a phase rj through the 
relation, cos 2r] = R + D cos cp, we obtain a solution for 
the quantity 7, 



= Sri - - + -KUs, s = ±, 



(A8) 



from which the energies of the Andreev bound states are 
found. 



,Acos(s?7 - 13/2), 



I, =sgn[sin(sr;-/?/2)]. 

(A9) 

The factor 6's = ±1 is fixed for each state by the condi- 
tion sin 7 > 0, which guarantees the exponential decay of 
the bound state wave functions into the superconducting 
leads. To simplify the further discussion, we assume the 
absence of spin selection, c?| = = r^. In this case, 

the relation D + R = 1 holds, and the parameter rj can 
be chosen as follows, 



/ — ■ V 
sin rj — V D sin — , 



(AlO) 



To write down an explicit form of the Andreev level 
wave functions, it is convenient to combine the envelopes, 
Eq. (jAl|l . in a four- vector, 0s — {4>t ^ then 



Fs 



Gs{x) 



(All) 
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,(x < 0) s6 



where 



and 



7 aJM->0), 



V2 



(A12) 



(A13) 



Fs = \ le- s^/D cos ^, 



Gsix) 



C ~T 1^ 



7^1-1 (A14) 



2hvFS 



^l- Dsm^{ip/2). 



APPENDIX B: SYMMETRY RELATIONS 

The symmetry relation, Eq. ([7|) generates the relations 
between the density matrices p and p, 

p(a;, x', t) ~ 5{x — x') — {iay)p*{x, x' , t){iay)\ (Bl) 

which extends to the respective single particle Hamilto- 
nians by virtue of the Liouville equation, 



h = —{iay)h*{iayy . 



(B2) 



Furthermore, Eq. (jB2p generates the symmetry relation 
between the respective eigen states. 



4>^{x) = {i(7y)(j)l{x), 



and eigen energies. 



(B3) 
iB4l 



Equations (|B3|) and (jB4[) establish mapping between the 
Hilbert spaces of the reciprocal Nambu representations, 
(^-representation and ^-representation. The eigen states 
of these representations form complete orthogonal sets, 
(01/, 01/') = {4>i^,4>i^i) = 5^v' ■ Furthermore, they obey an 
additional orthogonality relation, 



hy,4ly) = 0, 



(B5) 



that straightforwardly follows from the local identity, 
{(j}i,{x),iay(j)*{x)) = 0, the brackets here denote just a 
scalar product of two- vectors. 

The matrix elements of the reciprocal density matrices, 
p{x, x' , t) and x', t) in the respective eigen bases (cf. 
Eq. (HH), 



Puv'it) = {<j),y',p{t)<j)^), pu^j,{t) = {<t>^,,p{t)(f>i,), (B6) 
obey the symmetry relation, 

Puu' = <5^y' - pIi,' (B7) 
In particular, relation between the population numbers, 
fv = Puu and = pi,^, reads 



- 1 - 



(B8) 



The Andreev charge current and spin polarization are 
identical in both the representations, 

^■4 = ^ E ^^^^ (1 - 2/,) = ^ ^ d^E, (l - 2/,) , 



s=± 



Sa = ^(1 - /+ - /-) = -^(1 - /+ - /-). 



(B9) 
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